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Second Order Linear Differential Equation with 
Two Point Boundary Conditions (DOI Only) 

101 

Closed auxiliary for code Al| Code XA1, the floating 
point constant listing auxiliary, must also he used. 
The routine is entered as follows i 



p 


50 gF 
50 pF 


p + 1 


26 xF 



where g is the location of the auxiliary routine and 

x is the location of Code FA 1. 

See parameters . 

T *= (n + 1) (t + lKX>) + ^00 approximately where T is 

the total time in milliseconds > t is the time in 

milliseconds of the auxiliary routine and n is the 

number of sub -intervals into which the range of x 

is divided. 

9 decimal places or less, depending on the condition 

of the given equation. 

Preset parameters are S3, Sk, S5, S6, ST« 



CONTENT 
OOF 00 mF 

OOF OOJ0F 
OOF 00 kF 



USE 

m is the location of the floating point 
accumulator. 

J^ia the location of code Al 
Locations k + J (J = 0, ..., 2n +3) are 
used for temporary storage. The 
solutions y. ( i = 0, ... , n) will 
appear in the n + 1 successive locations 
k + a + 3 to k + 2n +:3- 



*»2» 



6 OQF 00 qf Locations q + j (J = 0, 1,...,10) con- 

tain the difference equation coefficients 
a i' ^i' c i* *** *** 8 differential equation 
coefficients j A,* B. * C , D. , ( i = 0, 
1, ... , n), the increment A (x) (8s6), 
and a temporary storage (9S6). The 
differential equation coefficients, A., 
B . , C. > D , calculated by the auxiliary 
ro routine, must he stored in locations 4s6, 

5S6, 6S6, and 7S6 respectively. 

7 OOF 00 ?F Locations r+J (j =0, 1, ..., 7) must 

contain the length of the range of x, L, 
and the boundary equation coefficients 
E, F, G, K, K, M expressed as floating 
point numbers, at the start of code 
FA 1- The quantities L, E, F, G, H, K 
and M are to be stored in locations S7* 
1S7, 2S7, 3S7> ^S7, 5S7, and 6S7 re- 
spectively. 
PROGRAM PARAMETER 

8 OQF 00 nF n is the number of sub*-intervals into 

•which the range of x is divided. 
DESCRIPTION Over an interval of length L at equally spaced points, 

x. (i = 0, 1, ..":••■« V&)s this routine will give the 
solution of the differential equation,, 







Ay* 


+ By 1 


+ cy 


= D, 


with 


the 


given boundary 


conditions , 






Sy 1 


+ Ty 


*■% 


*o 



and 

By* + Ky = M] 



x=x - nAx+x. 
m - 



•5- 



where A, B, C, D, E, F, 6, H, K and M are functions of x, the independent 
variable. 

An open subroutine, written using interpretive orders, 
must be provided starting at location g« It is to evaluate the coefficients 
A, B, C and D for a given x, x. < x < n Ax + x Q , of the differential 
equation* These calculated coefficients A, B, C and D must be stored 
successively starting at Location ks6* This auxiliary subroutine must be 
left vith the word 

8K IF 

83 (address of the 8th word of Code Fal) F. 
HOTE: Ax is available for use in the auxiliary subroutine. 

Ax is stored in 8s6. 

METHOD: This routine first calculates the coefficients of the 

difference equations, 



a i y i+ i + Yi + c i y i*i = d i 1 * °' n 



Vi + Vo _d o 

by + c y ,«d 
n rf n n n-1 n 

from the coefficients of the differential equations. The following 
equations are used to calculate a., b , c and d (0 < I < n): 

a = A + B^^ Ax/2 

b = C • Ax 2 - 2 A for i / 0, n 

c = A - B. Ax/2 Ax = L/n 

d, = D. Ax 2 
i 1 



-4- 



a o = A o E o 

b o - (E o c o - F o B o ) ' A x2/2 + A oV A x - A o E o 

c Q - 

d o J A o G o Ax+ ( E o D o" B o G o )Ax2 / 2 

a = 
n 

b = (KB - C H ) A x 2 /2 +AH+AK Ax 
n n n n n ' n n n n 

c = -A H 
n n n 

d = AM'Ax+(BM - H D ) A x 2 /2 
n n n n n n n ' 

At each step, these coefficients are stored in the four locations starting 
with location q. Then this routine evaluates the following formulae: 



h «■*■'■* 



i i "i-l 






■K-lrf: 



when i = 0, ... ,n and "V". = = "Z^ . p is calculated first and stored 

-1 -1 r 

in location k + 1. Then 1? is calculated and stored in location k + n + 3 
Finally ~Y~L is calculated and stored over JP in location k + 1. During 
the next time thru the loop^ 1 , 2^ and TT are calculated and stored in 
locations k + 2, k + n + 4 and k + 2 respectively. This process is con- 
tinued until , 2 s " and t~ have been calculated and stored. At this 
^n n n 

stage, the £1 are stored in locations k + 1 thru k + n + 1 and the "£7, 
are stored in locations k + n + 3 thru k + 2n + 3« When this has been 
done, the routine calculates the solution using the formula 



*i - ^i -*f2I + i i = n ' 



o. 



- 5 - 

where y =0. These solutions, y are stored over the tl in locations 
k + n + 3 thru k + 2n + 3. 

NOTE: Incorporated within the routine is a test to determine 

if £. = 0. When p. = 0, l/k x 10"^ is stored in place of this £ = 0. 

The method was adapted by S. Gill from a process used 
by L. Fox (L. Fox and H. H. Robertson; Automatic Digital Computation Symposium, 
National Physical Laboratory, England, March 1953. Paper 19, pp. 5, 6). 

Memo to FA 1 users. 

Code FA 1 is not compatable with SADOI. Consequently tape order should be as 
follows : 

1. DOI 

2. A 1 

3. XA 1 
k. FA 1 

5. Your program, written using DOI notation. (No symbolic addresses 
permis sable) . 

Important: Only five hole characters may be punched on the tape before DOI. 
The first non-5 hole character read must be the "8" in the first DOI instruction, 
Make sure that A 1 does not have a SADOI directive with a symbolic address at 
the beginning of it. (See A 1 writeup). 

RT: ] 2/9/59 

DATE February 16. 195^ 



-k_ 



PROGRAMMED BY R. Polivka 



APPROVED BY J. P. Nash 
ns 



LOCATION 


ORDER 


BOTES PAGE 1 FA 


'i 




1 < 
2 


*5 29F 
42 97L 
10 2QF 
42 7L 
36 83L 


" 

Set up link and address of 
auxiliary subroutine 






26 83L 








3 


4l S5 
50 3L 








4 


26 s4 
8K S8 








5 


8s 8s6 








6 


85 ST 

86 8s6 
8s 8s6 




-Calculate A x 




7 


8k IF 
83 (g)F 


From 33 j 
By 1 ^ 


to auxiliary subroutine 




8 

9 


/85 4s6\ 
^87 1S7/ 

8s s6 
81 5S6 


By 74, 78 '" 


replaced by 8K IF 

83 (54 )L 
Store a 




10 


87 2S7 
8s 9S6 








11 


85 6s6 








12 


87 1S7 
84 9S6 
87 8s6 








13 


87 8s6 
87 n+5+00 




Calculate the initial difference 




14 
15 


80 S6 
8s 9S6 
85 4s6 
87 8s6 




equation coefficients, a„, b , c 

* 0* 0* 
and d 




16 


87 2S7 
84 9S6 




1 





LOCATION 


ORDER 


NOTES PAGE 2 


. 17 


8S 1S6 
8K OF 






store b 


18 


8S 2S6 






store c ■ 





81 5S6 






19 


87 3S7 
8fe 986 








20 


85 1S7 
87 7S6 








21 


8k 9S6 
87 8s6 








22 


87 8s6 
87 N+5+00 






- 


23 


8s 9S6 
85 ks6 








2k 


87 8s6 
87 3S7 








25 


8k 9S6 
8s 3S6 






store d 



26 


6M jkL 
OK S8 






27 


81 2S6 
07 S5 


From 64,53 






28 


8k 1S6 
8j 69L 








29 


OS 1S5 






store w. 

Calculate p , C. , and *\/ . 


30 


81 2S6 

07 (k+n+2) 


F 






8k 3S6 








31 


06 1S5 
OS (k+n+3) 


F 




store 77 


32 


85 S6 

06 1S5 






X 



FA 1 



LOCATION 


ORDER 


NOTES PAGE 3 


53 


OS 1S5 
03 7L 






storet^ 


3* 


8j 77L 








83 (7)L 


by 82 






55 


/5k of\ 

18S S6# 






replaced by 8K IF 
82 6kL 


56 


81 ks6 
87 kS7 








37 


8s 2S6 
81 6s6 








38 


87 4S7 
8s 9S6 








39 


85 5S6 
87 5S7 








ko 


Qk 9S6 
87.836 








kl 


87 8s6 
87 N+5+00 








k2 


80 2S6 
8s 9S6 








k3 


85 ks6 
87 8s6 








kk 


87 5S7 
Qk 9S6 








^ 


8s 1S6 
81 4S7 






Calculate to the difference 

equation coefficients a , b 

n n 


k6 


87 7S6 
8s 9S6 






c and d . 
n n 


k-j 


85 5S6 
87 6S7 








kQ 


Qk 9S6 










87 8s6 









FA 1 



LOCATION 


ORDER 


NOTES PAGE k 


^9 


87 8S6 
87 N+5+00 










50 


8s 9S6 
85 ks6 










51 


87 8s6 
87 6S7 










52 


84 9S6 

85 3S6 










53 


8k IP 

83 27L 










5^ 


85 5S6 
87 8s6 


from 8 








55 


87 N+5+00 
8s 9S6 










56 


84 k$6 
8s S6 






store a 




57 


81 9S6 
84 ks6 










58 


8S 2S6 
85 6s6 






store 

calculate the difference 




59 


87 8S6 
87 8S6 






jequation coefficients a., 

c and d. i ^ 0,n. 
i 1 


V 


60 


80 4s6 
80 4s6 










61 


8s 1S6 
85 7S6 






store b. 

1 




62 


87 8s6 
87 8s6 




' 






63 


8s 3S6 
8k IF 






store d. 

1 




6* 


83 27L 






1 






IK s8 


From 35 








65 


11 (k+2n+3)F 
17 (k+n)F 1 









FA 1 



LOCATION 


ORDER 


NOTES PAGE 5 FA. 1 


66 


1^ (k+2n+2)F 




Calculate the solutions, y. 




IS (k+2n+2)F 




i 


67 


IL 1022F 
13 65L 






66 


8j 9^L 




' 




32 9^L 






69 


L3 S3 
32 70L 


from 28 




TO 


26 293V 








19 IF 




■ •Test for ft= 


71 


kO S3 






L5 73L 






72 


^0 1S3 
26 29S4 






73 


00 F 








00 28F 






7^ 


L5 76L 


Prom 26 






kO 8l 




Replace word 8 


75 


26 29S4 
00 7L 




by 8K IF 83 5^L 


76 


8k IF 








83 5^L 






77 


^9 S3 
L5 93L 


From 3^ 




78 


k2 8l 




Switch, 1st time go to part of 


79 


LI 98L 

ko 98L 




code to calculate a , b , c , d . 

n n n n 

2nd time, 




36 29Bh 




»- 


80 


L5 100L 




Go to part of code to calculate 




kO 35L 




the y 1 


81 


L5 93L 


' 






k2 3^L 







LOCATION 


ORDER 


: f 

NOTES PAGE 6 F/ 


82 


26 29S4 










00 F 






waste 


83 


L5 5F 
L4 8F 


From 2 






84 


42 65L 
F5 65L 








85 


L4 JOL 
42 92L 








86 


00 20F 
46 30L 






Set addresses in routine dependent 


87 


F5 65L 
F4 70L 






on the program parameters. 


88 


42 31L 
L5 92L 








89 


l4 8f 
42 66l 








90 


00 20F 
k6 66L 








91 


F5 66l 
00 20F 








92 


k6 65L 

kl (k+n+2)F 


5y 85 




* 


93 


26 3L 










oo 35L 






9* 


L3 98L 
40 8l 








95 


L5 99L 
40 35L 






Reset words 


96 


L5 75L 
42 34L 






-8, 34, and 35 


97 


32 97L 










22 ( )F 


1*0 


link 



LOCATION 



98 



99 



100 



■ njnin —»ww 



ORDER 



85 ^6 
87 1S7 
8k of 
8s s6 
8k if 
82 6kL 



NOTES 



PAGE 7 



FA 1 



wiwpa wi ap **mt M » i » w >ii t : 
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